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Abstract: Multijet plus missing energy searches provide universal coverage for theo- 
ries that have new colored particles that decay into a dark matter candidate and jets. 
These signals appear at the LHC further out on the missing energy tail than two-to- 
two scattering indicates. The simplicity of the searches at the LHC contrasts sharply 
with the Tevatron where more elaborate searches are necessary to separate signal from 
background. The searches presented in this article effectively distinguish signal from 
background for any theory where the LSP is a daughter or granddaughter of the pair- 
produced colored parent particle without ever having to consider missing energies less 
than 400 GeV. 
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1. Introduction 

A new energy frontier is being explored with the commissioning of the LHC, allowing 
for the best chance to discover physics beyond the Standard Model in over a decade. 
One of the most robust signals of new physics is anomalous contributions to jets plus 
missing energy from new strongly interacting particles that subsequently decay into 
jets and an invisible particle that may be the dark matter. 

In advance of the LHC, there have been many studies of the discovery potential in 
jets and missing energy [1-3]. Most of these studies study how to discovery a specific 
theory, e.g.the mSUGRA parameterization of the minimal supersymmetric Standard 
Model (MSSM) [4] (see [5-9] for other specific supersymmetric spectra), Universal Extra 
Dimensions (UED)s [10] or Little Higgs theories [11]. These searches present generally 
optimistic results for each specific theory, but different theories, even within the SSM 
may exhibit different signatures, [12]. How broadly applicable the searches optimized 
for specific theories are is an open question. This article begins to answer this question 
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by studying which searches are necessary to provide universal coverage for a broad class 
of spectra. 

Model-independent searches at hadron colliders are experimentally attainable by 
focusing on deviations from the Standard Model, and by parametrizing the details 
of a model in term of relevant observables like the masses of the particles produced 
[13]. [14,15] constructed a general search strategy at the Tevatron for discovering new 
physics in jets and missing energy ($r) regardless of the spectrum of new physics. The 
search strategy required combining several different search channels: lj + $t , 2j + $t , 
3j + $t and 4 + j + $t, where the last channel was inclusive. Inside each of these search 
channels, it was necessary to perform many measurements varying the visible energy 
(H T ) and missing energy. Several dozen measurements were necessary to fully cover 
the Tevatron's discovery potential. The primary challenge for the Tevatron is that new 
particles that the Tevatron can produce have similar masses to the Standard Model 
background processes. At the LHC, the new particles that can be searched for are 
typically much heavier than the SM background processes (dominantly top pair pro- 
duction) and therefore are more energetic. If there are new low mass colored particles 
that are similar in mass to the top, they are produced in such copious quantities that 
they are difficult to miss. This difference between the Tevatron and the LHC results 
in a much simpler search strategy and the possibility of a missed discovery is greatly 
reduced. 

The search strategy presented in this article will provide coverage for nearly all 
new colored particles that can be discovered in jets plus missing energy searches. This 
search strategy only uses the multijet (4 + j + $t) channel; the efficacy of this strategy 
originates from highly energetic initial state radiation (ISR) and final state radiation 
(FSR) jets that accompany new colored particle production at y/s = 14TeV at the 
LHC [16]. Parton shower/matrix element (PS/ME) [18] matching is used on both the 
background and signal in order to more accurately estimate the spectrum of additional 
hard radiation [19]. PS/ME matching increases the jet multiplicity of the signal, but it 
also increases the missing energy in signal events and therefore the discovery potential 
of the LHC. 

The efficacy of this search strategy is demonstrated through three benchmark su- 
persymmetric field contents: 

• Pair-produced gluinos in the supersymmetric Standard Model (SSM). Each gluino 
directly decays to two light-flavored quarks and the LSP. The LSP is the daughter 
of the original pair produced particle. 

• Pair-produced gluinos that cascade decay via a wino and jets with the wino 
subsequently decaying to a W ± and the LSP. The LSP is the granddaughter of 
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the original pair-produced particle. 

• Pair-produced squarks 1 in the SSM that decay to a quark and the LSP. 

These benchmark models are illustrative of the wider class of theories that produce jets 
plus missing energy signatures beyond supersymmetric and apply to UEDs and other 
theories which have new colored particles that decay into a dark matter candidate. 
The search strategies presented in this article are equally effective at separating gluinos 
that decay into multiple hard jets plus missing energy from background as they are at 
separating squarks that decay into a few soft jets plus missing energy from background. 

The organization of the article is as follows. Sec. 2 presents the general strategy 
for model independent jets plus missing energy searches. Sec. 3 describes the methods 
used to calculate the SM backgrounds and new physics signals. Sec. 4 shows the 
application of the search strategies to several minimal models. Finally, Sec. 5 discusses 
the types of physics that are not effectively searched for with the strategies described 
in this article. 

2. Search Strategy 

The goal of this article is to demonstrate that a simple search strategy will discover 
a broad class of new physics models without any loss of sensitivity to unexpected 
theories. In [15], the searches are first divided into mutually exclusive event samples 
with different numbers of jets: 

lj + $T 2]+$T 3]+$T 4 + J+$t (2.1) 

where 4 + j + fa is an inclusive, four or more jets plus missing energy search. At the 
Tevatron it is necessary to perform a series of searches in each of the multiplicities above 
in order to maximize the discovery potential for the full range of potential signals. 

This article will demonstrate that a series of measurements using only the 4 + j + fa 
channel is sufficient to gain broad sensitivity to essentially all theories that are visible in 
jets plus missing energy searches at the LHC. The effectiveness of the multijet search, 
even in theories with nig ~ m x o, ^ s a remarkable feature of a y/s = 14TeV LHC and 
requires that initial and final state radiation produce many, and occasionally all, of the 
jets used in the multijet search. 

The origin of this universality of the multijet search arises from several factors. The 
first is that radiation allows access to different components of the parton distribution 
functions (pdfs). Pair-produced new particles originate from the gg or qq initial states; 

1 The 9 flavors of squarks (9i,u£,d°) are degenerate. 
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however, an additional initial state jet allows access to the qg pdf that dominates over 
gg above v§ ~ 700 GeV. With two additional final state jets, pair production has 
access to the qq pdf and dominates above Vs ~ 2 TeV. Therefore, including additional 
jets in the signal allows access to more energetic pdfs which in turn allows for more 
energetic events. 

For compressed spectra, there is a kinematic reason why addition radiation is so 
important for distinguishing signal from background. Consider the case where a jet 
with energy Ej is radiated off a massive particle, X, of mass m. The four momenta of 
the jet and the massive particle are given as 

pj = (Ej, Ej) fx = (y/m* + E%, -E s ) . (2.2) 
The additional cost of energy for radiating a jet off a massive particle is given by 

Ei m > Ei 



5y/Ts s = Ej + Jw? + E]-m = \ 3 3 . (2.3) 



2Ej m < Ej 



Conversely, when radiating a jet off the SM backgrounds, letting m — > y/s in Eq. 2.3, 
with Ej ~ 2rrig 3> \fs B gives 



Sy/S B = Ej + yJsB + E] -VTs b = 2Ej. (2.4) 
Comparing signal and backgrounds yields 

yf§ s+j ~V§ s + Ej y/§ B+j ~2Ej + -^-; (2.5) 

letting v§s+j — V^B+j equalizes the CM energy for signal and background, thus 
increasing S/B. A similar procedure leads to the same relation in Eq. (2.5) when the 
jet comes from initial state radiation. The multitude of jets greatly simplifies the search 
strategy and reduces the number of measurements necessary to gain broad coverage of 
theories beyond the Standard Model. 

The remaining portion of this section is organized as follows. Sec. 2.1 classifies 
the search channels into different jet multiplicities and addresses the cuts on jet prs. 
Sec. 2.2 introduces the other kinematic and event shape variables that are useful at 
separating signal from background and specifies the series of searches that comprise 
the full search strategy. 

2.1 Classification of Search Channels 

In order to both trigger and ensure that the events are sufficiently distinctive compared 
to poorly understood backgrounds, the searches considered in this article have tighter 



-4- 





Selection Criteria for Different Jet Multiplicity Searches 




1 3 + $T 


2 j + fi T 


3 j + $ r 


4+j + $ T 


jl 


> 400 GeV 


> 1 00 GeV 


> f 00 GeV 


> 1 00 GeV 


h 


< 50 GeV 


> 50 GeV 


> 50 GeV 


> 50 GeV 


h 


< 50 GeV 


< 50 GeV 


> 50 GeV 


> 50 GeV 


J4 


< 50 GeV 


< 50 GeV 


< 50 GeV 


> 50 GeV 



Table 1: The selection criteria for the different jets plus missing energy samples. The 
abundance of ISR and FSR means that only the 4 + j + fix channel is necessary to maximize 
sensitivity. Jet pt cuts in GeV for the different exclusive search channels. The following cuts 
are applied: \rjj\ < 2.5, p T ,e < 20 GeV, A(j}(fi T ,ji) > 0.2 rad. (i = 1,2,3). 

cuts than those considered in the published search strategies from ATLAS [20] and 
CMS [21]. The cuts used in this article are loosely based on the benchmark ATLAS 
multijet + $t searches and the modifications to the searches tighten the searches [22]. 
This article finds that requiring a larger missing energy criteria is more effective at 
separating signal from background. 

The first step in devising a model- independent, inclusive search strategy is to clas- 
sify events into different jet multiplicities. The primary variables used are the pts of 
the jets. Each multiplicity has a definition of a high p?, central primary jet. After 
the primary jet there are secondary jets that can be central or forward with lower p? 
requirements. To avoid events falling between searches, a classification scheme that is 
relatively insensitive to the exact values used in the division of events is necessary. This 
is accomplished by selecting a common definition of a secondary jet and then classifying 
events by the number of secondary jets. The pt selection of the jets is as shown in 
Table I. 

There is a trade-off between having a higher versus lower requirement for jet pts. 
Higher jet pts mean that the jets are more likely to come from final state decay products 
and less likely to be produced in parton showering and radiation. Having a tighter 
requirement on jet Pts will reduce the SM backgrounds significantly. However, new 
physics may not have a widely spaced spectrum and the resulting final state jets may 
be soft. Therefore, by raising the requirements on jet pts, there is a loss of sensitivity to 
degenerate spectra. At the same time, the search strategy should pick out spectacularly 
hard jets from a sea of soft jets. In Sec. 2.2, other event variables, namely or M e g, 
are used to separate very energetic events. 

An alternative strategy would be to incorporate multiple searches using different 
requirements on jet pts in different searches. While this approach will increase reach, 
it does so at the expense of rapidly increasing the number of searches necessary to 
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cover the full range of possible signals. This alternative was studied in this article; 
ultimately, there is little gain in the discovery potential for this more complicated 
search. Therefore, keeping the jet pr requirements fixed and relatively modest is a good 
compromise and using Ht or M e s as a way of effectively increasing the Pt requirements 
on the jets. The value of the lower multiplicity searches is diminished by having modest 
Pt requirements. The reduced value of lower multiplicity searches has a side benefit: a 
smaller set of searches is necessary to cover the full range of new physics possibilities. 
Thus, a multijet search will be most effective even if the decay products of the parent 
particle are effectively invisible and all the jets must be generated through radiation. 

2.2 Observables 

The next step in the search strategy is to use kinematic or event shape variables to 
discriminate signal from background. Minimizing the number of searches performed 
and the number of event shape variables used is desirable for several reasons. 

First, the full scope of beyond the Standard Model signals frequently populates 
most of the spectrum of these kinematic variables. It is necessary to examine multiple 
cut values for each kinematic variable used in order to not eliminate signal. For instance, 
[14] showed that the high H T cuts used in some jets plus missing energy searches at the 
Tevatron could eliminate potentially discoverable signals. Using more variables results 
in an increase in the number of searches necessary to discover generic beyond the 
Standard Model signals. An additional reason to avoid too many kinematic variables 
in a search is that higher order corrections to the signal and background may alter 
the expected sensitivity, leading to greater theoretical uncertainty in the searches. The 
more kinematic variables used, the greater this problem is exacerbated. 

Three event shape variables are useful at separating signal from background and 
are studied below. The squark and gluino benchmark models are used to test the 
efficacy of cuts on different event shape variables that enhance visibility of new physics 
over background. A series of three searches are necessary to have broad sensitivity to 
beyond the Standard Model physics. 

Intuitively, $t is an important event shape variable to use in discriminating signal 
from background. Moreover, additional resolving power can frequently be attained by 
using a second kinematic variable such as M e g or Ht- In addition to the tried-and-true 
kinematic variables, the sensitivity in multijet and $t searches is explored by using 
the event shape variable, q;rts; first introduced in [23] as an alternative to using a fixed 
$t cut for removing the QCD backgrounds in dijet and missing energy searches. q;rts 
was initially defined as 




(2.6) 
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where 2 are the transverse momentum of the leading two jets in the event, rrij 1 j 2 
is the invariant mass between the leading and subleading jets, and 9jij2 is the angle 
between them. 

In QCD dijet events, where the jets are predominantly back-to-back, Eq. 2.6 gives 
ttRTS=0.5. Fluctuations of a jet's energy reduce ckrts an d typical QCD events have 
QJrts < 0.5. A SUSY signal can extend beyond «rts of 0.5 because $t is not being 
created due to fluctuations of the jet energy [23]. This suggested that a cut of «rts > 0.5 
could be used in lieu of a $t cut in dijet events. However, as discussed in Sec. 2.1, 
that radiation typically produces several additional jets in addition to the signal and 
reduces the efficacy of the dijet search. 

Initial studies for dijet and missing energy searches at the LHC indicate that «rts 
is an effective event shape variable for dijets [24,25]; however, exclusive dijet searches 
are never the most effective discovery channel for the classes of signals considered in 
this article even using orts- «rts can be generalized [26] to events with multijets by 
assigning each of the jets to one of two groups of jets, or superjets (Ji and J 2 ), so 
that the ratio Ptj 1 /ptj 2 is minimized, Qrts is then calculated with Eq. 2.6 using the 
4-momentum of superjets J\ and J 2 rather than the leading jets. The grouping of jets 
into superjets is illustrated in Fig. 1. 
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Figure 1: A sample three jet (ji, j 2 , J3,) plus fix event with thepx ratios and «rts computed 
for each of the three superjet groupings. The grouping that minimizes the pt imbalance of 
the two superjets is chosen and orts is computed from these two superjets. In this case, the 
first grouping has the lowest momentum imbalance. 

As an illustrative example, consider a 700 GeV gluino directly decaying to a 350 
GeV LSP. The left panel of Fig. 2 shows a orts distribution for a simulated QCD 



-7- 



sample and the electroweak SM backgrounds with no cut (top panel); its shape 
retains the feature orts > 0.50 as an alternative to a $t cut in dijets and $t searches. 
However, many new physics scenarios do not have a significant tail beyond orts — 0.50. 
Using orts in lieu of a $t requirement was not sufficient; however, using cvrts and a 
moderate $t could improve the searches over a $t cut alone. The right panel of Fig. 2 
shows the orts distributions with a $t cut of 400 GeV. The signal begins to deviate 
from backgrounds at moderate orts — 0.20 and by cutting at orts > 0.20, the signal 
is discoverable with C = 1 fb _1 . 
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Figure 2: orts distributions of combined SM backgrounds (black) and a signal spectrum 
with a 700 GeV gluino decaying to a 350 GeV LSP (blue). The left panels shows the orts 
distribution with no $t requirement. The right panel shows the distributions after a fix > 
400 GeV is imposed. Using orts > 0.20 in conjunction with fiq* in a set of searches allows 
for discovery with the criteria in Sec. 4. 



Ultimately three sets of searches are needed to separate signal from background. 
For the gluino pair-production case, additional coverage along the degeneracy region, 
m x o < rrig is gained by considering a set of multiple cuts on Ht while applying a 
high $t > 600 GeV cut and no orts cut (orts > 0). The "High $t Searches" also 
maximizes the reach of a search for pair-produced squarks when rriq > m x o. A second 
set of searches, the "ckrts Searches," in H T keeping $ T > 400 GeV but requiring 
«rts > 0.2 picks up those signal points of high mass gluinos with 2m x o < rrig, and 
the large squark mass low LSP mass region. An additional set of searches, the "Base 
Searches," in Hj- while keeping a medium > 400 GeV cut and no orts cut, is 
necessary to achieve discovery reach for those theories where the energy goes into visible 
states. For example, consider a theory where gluinos are pair-produced, but cascade 
decay via a Wino to the LSP and jets. The selection criteria for each of the three sets 
of searches is illustrated in Table 2. Lowering the $t requirement does not enhance 
the discovery potential due to the steeply rising background and by $t > 600 GeV, 
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the backgrounds are down to O(10 fb). The SM backgrounds $t spectrum falls faster 
than a typical signal point, thus lower cuts on $t only decrease S/B. These searches 
are not exclusive and cannot be combined together; however, each set is particularly 
suited for enhancing coverage in a region in the (rrig, m x o) mass parameter space. These 
strategies are applied to estimate the discovery sensitivity at the LHC at yfs = 14 TeV 
for pair-produced gluinos that either cascade or directly decay to the LSP and two jets, 
and to pair-produced squarks that decay directly to the LSP and a jet. 







«RTS 




Base 


> 400 GeV 


> 0.0 


> 600 GeV, 900 GeV, 1200 GeV, 1500 GeV 


High $ T 


> 600 GeV 


> 0.0 


> 600 GeV, 900 GeV, 1200 GeV, 1500 GeV 


«RTS 


> 400 GeV 


> 0.2 


> 600 GeV, 900 GeV, 1200 GeV, 1500 GeV 



Table 2: The different multijet plus missing energy searches used. The left table shows the 
division of searches into the "High fir" search, and "orts" search. The right table shows 
the different supplementary Ht cuts necessary to gain full coverage. The larger Ht cuts are 
effective at looking for highly non-degenerate spectra. 



3. Signal and Background Calculation 

This section describes the calculation of signal and background. This article uses 
leading order Monte Carlo calculations for processes up to 2 — > 4. The cross sections 
are scaled to the NLO inclusive cross section. 

The dominant backgrounds to $t + jets searches are tt + nj, W ± + nj, Z° + nj, 
and QCD. Subdominant contributions from W ± W ± , Z°W ± , single top, and Z°Z° are 
not considered in this article. 

Additional jets are calculated through two distinct approximation techniques in 
the Monte Carlo event generation. The first method uses Feynman diagram techniques 
with larger multiplicity final states. This approximation procedure gets the wide-angle, 
hard emission correct, but is computationally expensive and does not resum larger 
logarithms of kinematic variables. The parton shower is the second approximation 
procedure used to produce additional final state jets. Parton showers are appropriate in 
the collinear approximation and resum larger logarithms, but underestimate the amount 
of hard radiation. Both approximation schemes are used in Monte Carlo generators 
and matching final state jets to matrix element partons is necessary to avoid double 
counting a process represented in both methods. 
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The results of this study will rely on the 4 + j + $ T channel. Ideally one would 
want to generate up to four jets at matrix element level, but this is a computationally 
impractical task. Jets beyond those coming from the matrix element can be reasonably 
well approximated by the parton shower. For example, in [17], W + +jets were matched 
up to four jets. In this article, this background is matched up to three jets. The ratio 
of the inclusive W + + 4j rate to that of the inclusive W + + 3j rate from the samples 
generated for this study differ from the ratio from [17] by (9(15%). 

This article uses a MLM-based [18] k± shower matching scheme [19], where partons 
are clustered using the fcx-algorithm. The k± jet measure is defined as 

d 2 (t,j) = ARlmm( P 2 Ti .p 2 Tj ), (3.1) 

where 

AR^ = 2(cosh Arjij - cos A^-). (3.2) 

At the matrix element level, partons are required to be separated by a minimum cut-off 
qme [27] . After the parton shower and hadronization, the final state hadrons cluster 
into jets. Matching consists of identifying a jet with a parton. There are two instances 
to consider. First, for lower jet multiplicity samples, a jet is considered matched when 
it is within Q PS of a parton, namely rf(parton, j) < Q FS . If the jet cannot be matched, 
the event is discarded. Second, in the sample with the highest jet multiplicity (2 jets for 
ti and signal, 3 jets for W ± and Z°), additional unmatched jets are allowed provided 
that they are softer than the softest parton. The k± shower scheme sets Q ME = Q PS 
and in order to generate higher p? jets, it is necessary to set the minimum of the 
jets at parton level equal to Q ME . 

For the studies in this article, events were generated using Madgraph/MadEvent 
4.4.0 (4.4.2) [28] for each of the following background processes: 

pp W ± + nj 1 < n < 3 
pp^ Z° + nj 1 < n < 3 

pp^ti+nj < n < 2. (3.3) 

In order to have sufficient statistics at high missing energy, each of these backgrounds 
was calculated with different requirements of the pr of the heavy particle. For instance, 
the Z° + nj was calculated in four samples satisfying 



p TZ o < 100 GeV, 100 GeV < p TZ o < 200 GeV, 
200 GeV < p TZ o < 400 GeV, 400 GeV < p TZ o. (3.4) 



For each squark and gluino signal process considered, 100K events were produced: 

pp — >■ gg + nj < n < 2 

pp — >■ gg^ + nj < n < 2 (3-5) 

Pythia 6.4 is used to include initial and final state radiation beyond the radiation 
explicitly computed with Madgraph/MadEvent, and for the decays of all unstable par- 
ticles [29]. Finally, the ATLAS version of PGS 4 [30], with AR = 0.7, was used as a 
detector simulator. Appropriate K factors were included to ensure that the leading 
order in Madgraph/MadEvent production cross sections result in the correct total in- 
clusive production cross section. This article computed the K factors by comparing 
LO cross sections to theoretical predictions of the NLO cross section. As a cross check, 
after running the rescaled signal through PGS, different kinematic distributions that 
have been published by ATLAS were checked [22] . The calculations matched published 
distributions in [22] up to residual K factors that can be attributed to inaccurately 
modeled detector efficiencies such as lepton veto efficiencies in a multijet environment. 

As illustrated in Fig. 3, the generation of higher multiplicity jets from the parton 
shower underestimates the high p? and $t regions significantly, particularly in scenarios 
where rrig ~ m x o. In this region of parameter space, the jets coming from the decay of 
the gluino may be too soft to pass selection. Additionally, the two LSPs tend to line up 
back to back giving low $t- However, this topology is enhanced by including ISR jets. 
Emission of a hard ISR jet boosts the gluino pair and typically creates a momentum 
imbalance that gives high values of $t ■ 




1000 



#r (GeV) 

Figure 3: fir distributions for a matched (black) and unmatched (blue) 250 GeV gluino 
decaying to a 200 GeV LSP. 



- 11 - 



3.1 Background Subprocesses 

For both SM backgrounds and squark and gluino signals, CTEQ6L1 parton distribution 
functions [31] are used along with variable renormalization and factorization scales. 
This variable choice sets the scales to the central value of the transverse mass, my, 
of the event. At generation level at least one jet produced in association with X = 
tt, W ± , Z°, signal, must have pr > 100 GeV for each of the processes pp — > X + nj for 
n > 1. 

The backgrounds calculated in Sec. 3.1 are combined into a final estimate for the 
backgrounds in the three series of searches. The results are listed below in Table 3. 
This article uses a benchmark luminosity for the searches of 1 fb" 1 and therefore these 
searches should have relatively modest backgrounds. In Sec. 4, the sensitivity studies 
are performed and a systematic error of 30% is taken on all these cross sections. 

ti + nj 

ti+nj is the dominant background to three- and four-jets and fir searches at moderate 
values of missing energy (see Fig. 4). At matrix element level, the minimum transverse 
momentum of the jets is set to pr^ > 100 GeV = Q ME . Fig. 4 shows fir distributions 
for each SM background in the multijet channel. The LO cross section for tt + nj 
obtained is a = 483 pb. A K factor of 1.5 is taken from [32]. The ti M e ^ distribution 
is compared to that of [22], to obtain K e s,tt = 1.62. 




200 400 600 800 1000 

$ T (GeV) 

Figure 4: fix distributions for ti (blue), (green), Z° (orange), QCD (dashed red) and 
combined (black) in the multi-jet channel, after applying the selection criteria. 
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W ± + nj and Z° + nj 

W ± and Z° backgrounds are the most relevant backgrounds for high values of missing 
energy (see Fig. 4). At matrix element level, the minimum transverse momentum of the 
jets was set to prjet > 40 GeV, and sets the matching scale to Q ME = Q PS . Requiring 
that the px of the leading jet be at least 100 GeV does not conflict with the matching 
scheme if this cut is also placed at analysis level. To generate enough missing energy, 
only W and Z° decay modes with final state neutrinos are calculated: 

W ± ->e t + v t Z° -»■ vv. (3.6) 

Non-isolated muons are merged into the nearest jet and any objects PGS identifies as 
hadronic taus are converted into jets. 

The LO cross sections are 460 pb for W + + nj, 303 pb for W~ + nj, and 166 pb 
for Z° + nj. K factors are obtained from [33] for W ± and found to be 1.38 for W + 
and 1.35 for W~ . These results agree with those obtained by comparison with [34], 
and the effective K factor computed by matching the M e s to that from [22] . 

Z°+nj K factors are not available in the literature for sufficiently high multiplcities. 
In this article an effective K factor of 1.24 for Z° + nj is obtained by comparing the 
M eff distribution to that in [22]. 

QCD 

The QCD background gives dominant low missing energy contributions. The QCD 
background arises from a myriad of sources such as semileptonic decays of heavy quarks 
and mismeasurements or instrumental effects at detector level. Simulating QCD is 
challenging given the uncertainties from detector effects. A requirement on $ T > 
400 GeV is placed to remove the QCD backgrounds in the searches in this article. This 
cut is significantly stronger than many proposed searches that are frequently $t > 
200 GeV. The tighter requirement makes QCD sources of missing energy unimportant 
for the searches in this article. The QCD contributions for this article are produced in 
Madgraph/MadEvent 4.4.3; for the following parton level process: 

VV -> 3333 (3-7) 

Each jet is required to have > 50 GeV, and at least one of them must have 
Pt > 100 GeV. Additional jets come from showering and radiation done in Pythia. 
To generate enough statistics, the parton level events were divided into independent 
samples described by px cuts on the leading and sub-leading jets. The LO cross section 
obtained was oqcd = 72 nb. 
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The dominant contribution to $t from QCD arises from two sources: fluctuations 
of jet energies and neutrinos from heavy flavor decays. Jet energy fluctuations dominate 
low $t and can compete with ti, but is steeply falling. The QCD sample generated 
from this article used PGS to model the fluctuations jet energy to generate $t ■ Heavy 
flavor decays with real neutrinos dominate the tail of the QCD contribution to the $t 
distribution. No special attempt to quantify the heavy flavor contributions to missing 
energy; however, this tail takes over at I^t > 150 GeV, where QCD is subdominant 
to the high-pT sources of MET. All the searches in this article require $t > 400 GeV 
where the QCD/heavy flavor contributions to $t are well under control. 

4. Expected Sensitivity 



Standard Model Cross Section 


Ht 


Base 


High $t 


«RTS 


> 600 GeV 


155 fb 


25 fb 


63 fb 


> 900 GeV 


92 fb 


21 fb 


37fb 


> 1200 GeV 


40 fb 


12.6fb 


17fb 


> 1500 GeV 


16 fb 


6.2 fb 


8.0 fb 



Table 3: SM backgrounds cross sections (fb) in the multijet + $t channel performing the 
"High fi T " and the "orts searches". 

This section demonstrates the efficacy of the search strategy detailed in Table 2. 
The abundance of ISR and FSR jets means that one need only consider the 4 + j + for 
channel. The sensitivity to several minimal particle contents are shown in Sec. 4.1 to 
Sec. 4.3. This article performs counting experiments and discovery requires 

S > 5a tot (B) (4.1) 

where S is the number signal events, B is the number of background events, and a tot (B) 
is the combined statistical and systematic uncertainty in the background. 

The three classes of searches introduced in Table 2 are used to place discovery reach 
bounds for different sample theories. The first step is to calculate the cross section of 
the Standard Model backgrounds for each of the search channels. Table 3 shows the 
Standard Model predictions for the 3x4 searches. 

The uncertainty in any measurement is given by 

atot(B) = y / a Pois ( J B) 2 + (e Syst 5) 2 , (4.2) 
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Minimum Cross Section for 5 a Discovery 




Base 


High ft? 


«RTS 


> 600 GeV 


241 fb 


45 fb 


103 fb 


> 900 GeV 


146 fb 


39 fb 


63 fb 


> 1200 GeV 


68 fb 


25.9fb 


33 fb 


> 1500 GeV 


31 fb 


15.5 fb 


18.5 fb 



Table 4: The minimum cross section for a 5a discovery in the 4 + j + ft? channel when 
including Poisson fluctuations and a 30% systematic uncertainty. 

where B is the number of background events for a given measurement. The statistical 
error is Poisson with a variance, ap i S (B). ap Q i S (B) approaches \/B for large back- 
grounds. In addition to statistical fluctuations of backgrounds, there are systematic 
uncertainties that don't proportionally reduce with increased luminosity. es yst is the 
systematic uncertainty in the prediction and measurement and will ultimately limit the 
search reach. This systematic background represents both the theoretical uncertainty 
in the backgrounds and the experimental systematics. In the early running of the LHC 
a value of es ys t = 30% is a plausible estimate of the systematic uncertainty and is used 
throughout this article [35]. In order to discovery a signal with these requirements, a 
signal needs to have a cross section <j(S) > 5cr tot (B). These cross sections are listed in 
Table 4. 





Base 


High $t 


ttRTS 


> 600 GeV 


211 fb 


35 fb 


121 fb 


> 900 GeV 


137 fb 


30 fb 


72 fb 


> 1200 GeV 


62 fb 


19 fb 


31 fb 


> 1500 GeV 


26 fb 


10 fb 


13 fb 



Table 5: Cross sections in each search channel for a {nig, m x o) = (800 GeV, 300 GeV) sample 
spectrum. The Ht = 600 GeV, 900 GeV orts searches have sufficiently large cross sections 
to be discovered. 

Up to this point, the procedure is entirely model-independent. The next step is 
to take the model of choice and compute the number of events that it predicts for 
each of the searches described in Table 2. As an illustration of how a search for a 
signal point should proceed, Table 5 displays the multijet cross sections for a 800 GeV 
gluino decaying to a 300 GeV LSP. The greatest significance is found by applying a cut 
(H T > 600 GeV, orts > 0.2), the ojrts search. 

The remaining portion of this article demonstrates how these 3x4 searches listed 
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in Table 2 can effectively cover the kinematic parameter space of both jet-rich signals 
like gluino pair production and jet-poor signals like squark pair production. 

4.1 Gluino Sensitivity 

The model-independent search strategy described above is applied to a search for pair- 
produced gluinos. Each gluino decays via: 

g^q + q + X°- (4-3) 

The masses of the gluino and the LSP, the most relevant experimental observables, are 
unconstrained. The minimum transverse momentum of the ISR jets is set to px > 100 
GeV at parton level. Each gluino can decay to any one of the lightest four standard 
model quarks and its antiquark. The gluinos are produced dominantly through QCD 
interactions. The NLO cross sections were calculated in Prospino 2 . 1 [36,37] and used 
to obtain K factors for each gluino mass point generated. MadGraph is used to generate 
additional hard radiation which is important for degenerate spectra. The matching of 
matrix element radiation to parton shower jets is performed in an identical procedure 
to the SM tt backgrounds. 

The (mg,m x o) mass parameter space is partitioned into three interesting regions, 
and although these subsets of mass parameter space are not exclusive, they serve to 
illustrate the efficacy of the sets of searches listed in Table 2. First, the bulk region, 
where the mass splitting between parent and daughter particle is larger than the mass 
of the daughter particle and is a universal feature in mSUGRA models. In these 
benchmark spectra, the visible and missing energy of the jets coming from a SUSY 
decay are large enough for f/r.signai to beat the SM backgrounds. Thus this region is 
accessible to multijet and $t searches through both the high $t and the «rts searches. 

The second region of interest is the degenerate spectra scenario where the mass 
splitting between the parent and daughter particle is less than the daughter particle 
mass. Degenerate spectra can only be discovered because there are energetic ISR and 
FSR jets associated with their production. Radiation from heavy objects is suppressed, 
and therefore most of the jets in heavier degenerate spectra come from ISR and are at 
larger rapidity. The searches described in this article use a rapidity cut of \r)j\ < 2.5, 
which allows even forward jets to contribute to the signal. 

The degenerate gluino and LSP scenario has a slower falling distribution, as ISR 
and FSR jets populate this topology typically in regions of high 77. In non-degenerate 
spectra, as in the 800 GeV gluino, the more central jets, both from FSR and decay, play 
a greater role. Having hard radiation jets recoil off of a heavy gluino particle boosts the 
gluino pair and creates a momentum imbalance which typically gives large This 
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second scenario is discoverable through the high ft? search. Additional discovery reach 
is possible due to the large gluino production cross section and degenerate spectra with 
masses up to 500 GeV, which are visible at the LHC with this set of searches. 

The third subset of mass parameter of space to consider is when m x o is roughly 
within an additive factor of m g {e.g. m x o < 2m§). This is a region where the reach can 
be extended performing an ckrts search. The low ojrts region, qrts < 0.2, vetoed in the 
search is largely populated by SM backgrounds. Compressed spectra signals populate 
larger ojrts due to a greater energy mismatch in the jets from SM, pushing «rts from 
Eq. 2.6 to lower values. With signals however, the energetic jets from radiation are 
leveled in energy with the energy from the decay jets. For instance, a 750 GeV gluino 
decaying to a 350 LSP gives roughly 200 GeV jets. This pushes orts in Eq- 2.6 to higher 
values than 0.2. The results are summarized in Fig. 5, which shows the discovery reach 
at the LHC with yfs = 14TeV and an integrated luminosity of 1 fb _1 for pair-produced 
gluinos. 



600 




Gluino Mass (GeV) 

Figure 5: Discovery reach at the LHC for directly decaying gluino pair production at C = 
1 fb _1 . In light blue, the reach of a high Pt set of searches, and the additional reach due to 
an ckrts search (red). Also shown is the 2a exclusion at the Tevatron with 4 fb _1 of data 
(dark blue) from [15]. 

4.2 Cascade Decays 

The case studied so far assumes that each gluino decays directly to two jets and the LSP. 
It is important to consider an alternate scenario where the energy goes into intermediate 
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states, i.e. cascade decays. Here the gluinos are allowed to decay via: 

9 -> q + q + W W -> W* + X °- (4-4) 

In the mass spectrum, the Wino is assumed to lie halfway between the gluino and LSP, 
namely = 1/2 (m,g + rriLSp)- This spacing in the mass spectrum was shown in [14] 
to be the most pessimistic scenario at the Tevatron. Cascade decays lower the $t 
of the event because the energy is distributed to more visible states. These theories 
require lowering the $t to have maximal sensitivity. Fig. 6 shows the coverage by the 
three sets of searches introduced. The reach of the high $t search turns off when there 
is enough phase space for the W boson to go on-shell. 



500 




200 400 600 800 

Gluino Mass (GeV) 



Figure 6: Discovery reach at the LHC for cascade-decaying gluino pair production at C = 
1 fb^ 1 . In light blue, the reach of a high fix > 600 GeV set of searches, and the reach due to 
an orts search (red). The reach of a Base Search ft? > 400 GeV is shown in green. 



4.3 Squark Sensitivity 

The third class of theories studied in this article include spectra where the squarks 
are significantly lighter than the gluinos. In this case, pair-produced squarks are the 
dominant mechanism for discovering new physics. Squarks are scalars that are triplets 
of color, therefore, their production cross section and the associated radiation is smaller 
than gluinos. As a way of studying the efficacy of the multijet search strategies proposed 
in this article, the squarks are forced to decay directly to the LSP: 

q^q + X °. (4.5) 
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This topology indicates that a dijet plus missing energy search should be effective. 
However, ISR and FSR jets play a crucial role in this well. The matching scale 

was set to Q ME = Q PS = pxj > 75 GeV. A if-factor is used by calculating the NLO 
cross section in Prospino 2. 1 as in the case of gluinos [36,37]. Parton shower/matrix 
element matching is performed upon the signal to generate additional radiation. Due 
to the smaller color factors, there is significantly less overall radiation. Nevertheless, 
the additional jets in the event and the moderate pt requirement on them are sufficient 
to supplement the two final state jets from squark pair production. 

A squark signal point and the background $t distributions for squarks are shown 
in Fig. 7 for the multijet and dijet plus missing energy searches. While the overall 
event rate of dijet and $r is larger than the multijet, the discovery potential is greater 
for the multijet search. 




#r (GeV) $ T (GeV) 



Figure 7: fix distributions for SM backgrounds (black) and a 250 GeV q (blue) decaying 
to a 50 GeV LSP and a light flavored quark in dijet (left panel) and multijet (right panel) 
channels. 



This study demonstrates that the exclusively multijet search strategy presented 
in Sec. 2 and Sec. 2.2 is broadly applicable to many theories that do not have event 
topologies that appear dominantly as multijet events. As a comparison between the 
reaches of the multijet and dijet channels, Fig. 8 shows discovery reach plots of both 
multijet and dijet and missing energy searches for pair-produced squarks decaying to 
a jet and the LSP, with yfs = 14TeV and C = lfb - . Also shown is the reach of an 
inclusive dijet and $t search, where the p? of additional jets beyond the leading two 
is not vetoed. Moreover, increasing the pt requirement of the leading two jets in the 
dijet and $t search was also studied but found to be less sensitive than the criteria 
in Table 1, particularly in the degenerate spectra region. Scalar color-triplets have 
a smaller production cross section than fermionic color-octets and results in lowering 
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discovery potential at the LHC and is particularly difficult if the systematic errors are 
not reduced from 30%. 




4.4 Associated Squark Gluino Sensitivity 

As a final check, three benchmark points of squark-gluino associated production were 
computed: 

(m- g ,mq,m x0 ) = (400 GeV, 400 GeV, 200 GeV), 
(600 GeV, 400 GeV, 200 GeV), 
(250 GeV, 400 GeV, 150 GeV) . (4.6) 

In the first benchmark point, the gluino directly decays to two jets and the LSP, whereas 
the squark decays to a jet and the LSP. In the second case, the gluino decays to the 
squark and a SM quark; the squark, as in the previous case decays to a jet and the 
LSP. In the final benchmark, point the squark decays to the gluino and a quark; the 
gluino, in turn, decays to two jets and the LSP. Naively, the first and third benchmark 
points should be covered by trijet and missing energy searches, whereas the second 
point by a dijet and missing energy search. However, each of these benchmark spectra 
is separated from background by the multijet and missing energy searches depicted in 
Table 2. 
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5. Discussion 



Jets and missing energy searches are one of the most promising windows into new 
physics at the electroweak scale. The production of new colored particles comes with 
abundant energetic ISR and FSR jets. This article built upon this observation to 
propose a model-independent search strategy where nearly all classes of theories of 
beyond the SM containing colored particles which decay to jets and a neutral stable 
particle can be discovered in the 4 + j + $t channel at yfs = 14 TeV. The approach is to 
perform a set of three searches that vary the amount of visible energy, missing energy, 
or a RTS . 

This search strategy was applied to several classes of theories. The first considered 
pair-produced gluinos, with each gluino decaying to two jets and the LSP. The sec- 
ond class of theories included pair produced gluinos that decayed in a cascade through 
an intermediate particle and finally to the LSP. The third class of theories were pair- 
produced squarks where each squark directly decayed to a jet and the LSP. As a check, 
three benchmark spectra in the interpolating region of associated squark-gluino pro- 
duction were studied. 

The primary finding of this article is that incorporating parton shower matrix 
element matching to the signal increases the visibility of the signals by raising the 
missing energy in the event, particularly for heavy LSP masses. This results in higher 
missing energy cuts (greater than 400 GeV) being more effective at discovering new 
physics. 

Lower missing energy searches are not necessary until long cascades become normal. 
If the LSP is a granddaughter of the produced particle, a 400 GeV missing energy cut is 
most effective at discovering the spectrum. When the cascades become longer, so that 
the LSP is a great-granddaughter of the produced particle, the missing energy can be 
sufficiently degraded so that it is necessary to use lower missing energy cuts. A typical 
3-step cascade that has the LSP as the great granddaughter can still be discovered 
effectively with high missing energy cuts, but those where there is a hierarchal spectrum 
of masses such as 

rrig : : : rrig = 8:4:2:1 (5.1) 

where S is a singlino that appears in many NMSSM models. The difficulty in using 
high $t requirements arises because the momentum of the final LSP carries so little 
of the original particle's rest mass energy or original momentum. Spectra such as 

rrig : : : rrig = 8:7:6:1 or m g : : : m s = 16 : 4 : 3 : 2 (5.2) 

produce LSPs that still carry a reasonable fraction of the produced particle's momen- 
tum. In the first case, most of the center of mass energy is distributed in the last decay. 
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In the second case, most is distributed in the first decay and then because the sub- 
sequent decays each have daughter particles that are non-relativistic, the momentum 
from the original daughter particle is carried transferred to the granddaughter and then 
to the great granddaughter efficiently. 

Acknowledgements 

We would like to thank Johan Alwall, Mariangela Lisanti, and Jesse Thaler for helpful 
discussions. We thank Daniele Alves, Jared Kaplan, and Mariangela Lisanti for reading 
early versions of the draft and providing incredibly valuable feedback. EI and JGW 
thank the GGI institute for its hospitality in the later stages of this work. JGW is 
supported by the DOE under contract DE-AC03-76SF00515 and by the DOE's Out- 
standing Junior Investigator Award. 

References 

[1] J. Alitti et al. [UA2 Collaboration], Phys. Lett. B 235, 363 (1990). X. Portell [CDF Run 
II Collaboration], AIP Conf. Proc. 842, 640 (2006) [arXiv:hep-ex/0609017]. S. M. Wang 
[CDF and DO Collaborations], AIP Conf. Proc. 1078, 259 (2009). A. Abulencia et al. 
[CDF Collaboration], Phys. Rev. Lett. 97, 171802 (2006) [arXiv:hep-ex/0605101]; The 
DJ3 Collaboration, DJ3 CONF 4400 (2007). V. M. Abazov et al. [DO Collaboration], 
arXiv:0712.3805 [hep-ex]. S. M. Wang [CDF and DO Collaborations], AIP Conf. Proc. 
1078, 259 (2009). 

[2] S. Mrenna, G. L. Kane, G. D. Kribs and J. D. Wells, Phys. Rev. D 53, 1168 (1996) 
[arXiv:hep-ph/9505245]. H. Baer, C. h. Chen, F. Paige and X. Tata, Phys. Rev. D 52, 
2746 (1995) [arXiv:hep-ph/9503271]. W. Beenakker, R. Hopker, M. Spira and 
P. M. Zerwas, Nucl. Phys. B 492, 51 (1997) [arXiv:hep-ph/9610490]. I. Hinchliffe, 
F. E. Paige, M. D. Shapiro, J. Soderqvist and W. Yao, Phys. Rev. D 55, 5520 (1997) 
[arXiv:hep-ph/9610544]. B. K. Gjelsten, D. J. Miller and P. Osland, JHEP 0412, 003 
(2004) [arXiv:hep-ph/0410303]. B. K. Gjelsten, D. J. Miller and P. Osland, JHEP 0506, 
015 (2005) [arXiv:hep-ph/0501033]. K. Kawagoe and M. M. Nojiri, Phys. Rev. D 74, 
115011 (2006) [arXiv:hep-ph/0606104]. S. G. Kim, N. Maekawa, K. I. Nagao, 
M. M. Nojiri and K. Sakurai, JHEP 0910, 005 (2009) [arXiv:0907.4234 [hep-ph]]. 
M. M. Nojiri, K. Sakurai, Y. Shimizu and M. Takeuchi, arXiv:0808.1094 [hep-ph]. 
D. E. Morrissey, T. Plehn and T. M. P. Tait, arXiv:0912.3259 [hep-ph]. P. Nath et al, 
arXiv: 1001. 2693 [hep-ph]. 

[3] C. Macesanu, C. D. McMullen and S. Nandi, Phys. Rev. D 66, 015009 (2002) 

[arXiv:hep-ph/0201300]. H. C. Cheng, K. T. Matchev and M. Schmaltz, Phys. Rev. D 



- 22 - 



66, 056006 (2002) [arXiv:hep-ph/0205314]. C. Macesanu, C. D. McMullen and S. Nandi, 
Phys. Lett. B 546, 253 (2002) [arXiv:hep-ph/0207269]. 

[4] S. Dimopoulos and H. Georgi, Nucl. Phys. B 193, 150 (1981). H. P. Nilles, Phys. Rep. 
110, 1(1984). B. K. Gjelsten, D. J. Miller and P. Osland, JHEP 0412, 003 (2004) 
[arXiv:hep-ph/0410303]. B. K. Gjelsten, D. J. Miller and P. Osland, JHEP 0506, 015 
(2005) [arXiv:hep-ph/0501033]. 

[5] S. Dimopoulos, S. D. Thomas and J. D. Wells, Nucl. Phys. B 488, 39 (1997) 

[arXiv:hep-ph/9609434]. G. F. Giudice and R. Rattazzi, Phys. Rept. 322, 419 (1999) 
[arXiv:hep-ph/9801271]. 

[6] K. Choi, A. Falkowski, H. P. Nilles, M. Olechowski and S. Pokorski, JHEP 0411, 076 
(2004) [arXiv:hep-th/0411066]; K. Choi, A. Falkowski, H. P. Nilles and M. Olechowski, 
[arXiv:hep-th/0503216]. W. S. Cho, Y. G. Kim, K. Y. Lee, C. B. Park and Y. Shimizu, 
JHEP 0704, 054 (2007) [arXiv:hep-ph/0703163]. K. Choi, A. Falkowski, H. P. Nilles and 
M. Olechowski, Nucl. Phys. B 718, 113 (2005) [arXiv:hep-th/0503216]. A. Falkowski, 
O. Lebedev and Y. Mambrini, JHEP 0511, 034 (2005) [arXiv:hep-ph/0507110]. H. Baer, 
E. K. Park, X. Tata and T. T. Wang, JHEP 0608, 041 (2006) [arXiv:hep-ph/0604253]. 
H. Baer, E. K. Park, X. Tata and T. T. Wang, Phys. Lett. B 641, 447 (2006) 
[arXiv:hep-ph/0607085]. H. Baer, E. K. Park, X. Tata and T. T. Wang, JHEP 0706, 
033 (2007) [arXiv:hep-ph/0703024]. W. S. Cho, Y. G. Kim, K. Y. Lee, C. B. Park and 
Y. Shimizu, JHEP 0704, 054 (2007) [arXiv:hep-ph/0703163]. 

[7] L. Randall and R. Sundrum, Nucl. Phys. B 557, 79 (1999) [arXiv:hep-th/9810155]; 

G. F. Giudice, M. A. Luty, H. Murayama and R. Rattazzi, JHEP 9812, 027 (1998) 
[arXiv:hep-ph/9810442]. T. Gherghetta, G. F. Giudice and J. D. Wells, Nucl. Phys. B 
559, 27 (1999) [arXiv:hep-ph/9904378]. 

[8] R. Kitano and Y. Nomura, Phys. Lett. B 631, 58 (2005) [arXiv:hep-ph/0509039]. 
D. Feldman, Z. Liu and P. Nath, Phys. Rev. D 80, 015007 (2009) [arXiv:0905.1148 
[hep-ph]]. T. Plehn and T. M. P. Tait, J. Phys. G 36, 075001 (2009) [arXiv:0810.3919 
[hep-ph]]. R. Kitano and Y. Nomura, Phys. Lett. B 631, 58 (2005) 
[arXiv:hep-ph/0509039]. H. Baer, S. de Alwis, K. Givens, S. Rajagopalan and 

H. Summy, arXiv:1002.4633 [hep-ph]. 

[9] C. H. Chen, M. Drees and J. F. Gunion, Phys. Rev. D 55, 330 (1997) [Erratum-ibid. D 
60, 039901 (1999)] [arXiv:hep-ph/9607421]. G. Anderson, H. Baer, C. h. Chen and 
X. Tata, Phys. Rev. D 61, 095005 (2000) [arXiv:hep-ph/9903370]. G. A. Blair, 
W. Porod and P. M. Zerwas, Eur. Phys. J. C 27, 263 (2003) [arXiv:hep-ph/0210058]. 

[10] T. Appelquist, H. C. Cheng and B. A. Dobrescu, Phys. Rev. D 64, 035002 (2001) 
[arXiv:hep-ph/0012100]. H. C. Cheng, K. T. Matchev and M. Schmaltz, Phys. Rev. D 



- 23 - 



66, 036005 (2002) [arXiv:hep-ph/0204342]. T. Appelquist, H. C. Cheng and 
B. A. Dobrescu, Phys. Rev. D 64, 035002 (2001) [arXiv:hep-ph/0012100]. 

[11] N. Arkani-Hamed, A. G. Cohen and H. Georgi, Phys. Lett. B 513, 232 (2001) 

[arXiv:hep-ph/0105239]; N. Arkani-Hamed, A. G. Cohen, T. Gregoire and J. G. Wacker, 
JHEP 0208, 020 (2002) [arXiv:hep-ph/0202089]. N. Arkani-Hamed, A. G. Cohen, 
E. Katz, A. E. Nelson, T. Gregoire and J. G. Wacker, JHEP 0208, 021 (2002) 
[arXiv:hep-ph/0206020]; N. Arkani-Hamed, A. G. Cohen, E. Katz and A. E. Nelson, 
JHEP 0207, 034 (2002) [arXiv:hep-ph/0206021]. H. C. Cheng and I. Low, JHEP 0309, 
051 (2003) [arXiv:hep-ph/0308199]. H. C. Cheng and I. Low, JHEP 0408, 061 (2004) 
[arXiv:hep-ph/0405243]. E. Katz, J. y. Lee, A. E. Nelson and D. G. E. Walker, JHEP 
0510, 088 (2005) [arXiv:hep-ph/0312287]. H. C. Cheng, I. Low and L. T. Wang, Phys. 
Rev. D 74, 055001 (2006) [arXiv:hep-ph/0510225]. A. Freitas and D. Wyler, JHEP 
0611, 061 (2006) [arXiv:hep-ph/0609103]. T. Gregoire and E. Katz, JHEP 0812, 084 

(2008) [arXiv:0801.4799 [hep-ph]]. 

[12] C. F. Berger, J. S. Gainer, J. L. Hewett and T. G. Rizzo, JHEP 0902, 023 (2009) 

[13] N. Arkani-Hamed, P. Schuster, N. Toro, J. Thaler, L. T. Wang, B. Knuteson and 
S. Mrenna, arXiv:hep-ph/0703088. T. Aaltonen et al. [CDF Collaboration], 
arXiv:0712.2534 [hep-ex]. T. Aaltonen et al. [CDF Collaboration], arXiv:0712.1311 
[hep-ex]. C. Henderson [CDF Collaboration], arXiv:0805.0742 [hep-ex]; G. Choudalakis, 
PhD thesis, Massachusetts Institute of Technology (2008); N. Arkani-Hamed, 
G. L. Kane, J. Thaler and L. T. Wang, JHEP 0608, 070 (2006) [arXiv:hep-ph/0512190]. 

[14] J. Alwall, M. P. Le, M. Lisanti and J. G. Wacker, Phys. Rev. Lett. B 666, 34 (2008) 
[arXiv:0803.0019 [hep-ph]]. 

[15] J. Alwall, M. P. Le, M. Lisanti and J. G. Wacker, arXiv:0809.3264 [hep-ph]. 

[16] J. Alwall, K. Hiramatsu, M. M. Nojiri and Y. Shimizu, Phys. Rev. Lett. 103, 151802 

(2009) [arXiv:0905.1201 [hep-ph]]. 

[17] J. Alwall et a/., Eur. Phys. J. C 53, 473 (2008) [arXiv:0706.2569 [hep-ph]]. 

[18] M. H. Seymour, Comput. Phys. Commun. 90, 95 (1995) [arXiv:hep-ph/9410414]. 
S. Mrenna and P. Richardson, JHEP 0405, 040 (2004) [arXiv:hep-ph/0312274]. 
M. L. Mangano, M. Moretti, F. Piccinini, R. Pittau and A. D. Polosa, JHEP 0307, 001 
(2003) [arXiv:hep-ph/0206293]. A. Schalicke and F. Krauss, JHEP 0507, 018 (2005) 
[arXiv:hep-ph/0503281]. F. Krauss, A. Schalicke, S. Schumann and G. Soff, Phys. Rev. 
D 72, 054017 (2005) [arXiv:hep-ph/0503280]. 

[19] T. Plehn, D. Rainwater and P. Z. Skands, Phys. Lett. B 645, 217 (2007) 

[arXiv:hep-ph/0510144]. J. Alwall, S. de Visscher and F. Maltoni, JHEP 0902, 017 



- 24 - 



(2009) [arXiv:0810.5350 [hep-ph]]. S. de Visscher, J. Alwall and F. Maltoni, AIP Conf. 
Proc. 1078, 293 (2009). J. Alwall, K. Hiramatsu, M. M. Nojiri and Y. Shimizu, Phys. 
Rev. Lett. 103, 151802 (2009) [arXiv:0905.1201 [hep-ph]]. 

[20] J. Abdallah et al. [ATLAS Collaboration], ATL-COM-PHYS-2009-189. N. Ozturk 
[ATLAS Collaboration and CMS Collaboration], arXiv:0910.2964 [hep-ex]. 
A. W. Phillips, CERN-THESIS-2008-105. T. Eifert [ATLAS Collaboration], 
ATL-COM-PHYS-2008-193. D. Costanzo, J. Phys. Conf. Ser. 110, 072007 (2008). 
M. Uslenghi, CERN-THESIS-2008-020. 

[21] G. Lungu [CMS Collaboration], arXiv:0910.3310 [hep-ex]. N. Ozturk [ATLAS 
Collaboration and CMS Collaboration], arXiv:0910.2964 [hep-ex]. C. Hof [CMS 
Collaboration], CERN-CMS-CR-2008-106. A. Tricomi [CMS Collaboration], J. Phys. 
Conf. Ser. 110, 062026 (2008). T. Yetkin and M. Spiropulu [CMS Collaboration], Acta 
Phys. Polon. B 38, 661 (2007). 

[22] G. Akimoto et al. [ATLAS Collaboration], ATL-COM-PHYS-2009-187. 

[23] L. Randall and D. Tucker-Smith, Phys. Rev. Lett. 101, 221803 (2008) [arXiv:0806.1049 
[hep-ph]]. 

[24] H. Baer, V. Barger, A. Lessa and X. Tata, JHEP 0909, 063 (2009) [arXiv:0907.1922 
[hep-ph]]. 

[25] CMS Collaboration, CMS PAS SUS-08-005 (2008). U. De Sanctis [ATLAS 
Collaboration and CMS Collaboration], Nuovo Cim. 123B, 793 (2008). 

[26] Private communication with L. Randall (June 8, 2009) at SUSY'09, Boston, MA, USA. 

[27] S. Catani, Y. L. Dokshitzer, M. H. Seymour and B. R. Webber, Nucl. Phys. B 406, 187 
(1993). 

[28] J. Alwall et al, JHEP 0709, 028 (2007) [arXiv:0706.2334 [hep-ph]]. 

[29] T. Sjostrand, S. Mrenna and P. Skands, JHEP 0605, 026 (2006) 
[arXiv:hep-ph/0603175]. 

[30] J. Conway, PGS: Pretty Good Simulator, http://www.physics.ucdavis.edu/ conway/ 
research / software /pgs/ pgs4-gener al . htm . 

[31] J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadolsky and W. K. Tung, JHEP 
0207, 012 (2002) [arXiv:hep-ph/0201195]; D. Stump, J. Huston, J. Pumplin, 
W. K. Tung, H. L. Lai, S. Kuhlmann and J. F. Owens, JHEP 0310, 046 (2003) 
[arXiv:hep-ph/0303013]. 

[32] N. Kidonakis and R. Vogt, Phys. Rev. D 78, 074005 (2008) [arXiv:0805.3844 [hep-ph]]. 



- 25 - 



[33] C.F. Berger, Z. Bern, L.J. Dixon, F. Febres Cordero, D. Forde, T. Gleisberg, H. Ita, 
D.A. Kosower, D. Maitre, [arXiv:0907.1984 [hep-ph]]; C.F. Berger, Z. Bern, L.J. Dixon, 
F. Febres Cordero, D. Forde, T. Gleisberg, H. Ita, D.A. Kosower, D. Maitre, 
[arXiv:0909.4949 [hep-ph]]. 

[34] J. Campbell, R.K. Ellis, D. Rainwater, Phys. Rev. D 68, 094021 (2003) 
[arXiv:hep-ph/0308195]. 

[35] G. Polesello, (Oct. 21, 09). "Susy at ATLAS with Early LHC Data." Seminar at 

Searching for New Physics at the LHC, Galileo Galilei Institute for Theoretical Physics, 
Arcetri, Italy. e syst = 50% stated in slides, e syst = 30% stated at presentation. 

[36] W. Beenakker, R. Hopker and M. Spira, arXiv:hep-ph/9611232. 

[37] E. Mirabella, JHEP 0912, 012 (2009) [arXiv:0908.3318 [hep-ph]]. A. Kulesza and 
L. Motyka, Phys. Rev. D 80, 095004 (2009) [arXiv:0905.4749 [hep-ph]]. U. Langenfeld, 
AIP Conf. Proc. 1200, 333 (2010) [arXiv:0908.3567 [hep-ph]]. 



- 26 - 



